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ABSTRACT 

Non-thermal motions in the intra-cluster medium (ICM) are believed to play a non- 
negligible role in the pressure support to the total gravitating mass of galaxy clusters. 
Future X-ray missions, such as ASTRO-H and ATHENA, will eventually allow us 
to directly detect the signature of these motions from high-resolution spectra of the 
ICM. In this paper, we present a study on a set of clusters extracted from a cosmolog- 
ical hydrodynamical simulation, devoted to explore the role of non-thermal velocity 
amplitude in characterising the cluster state and the relation between observed X-ray 
properties. In order to reach this goal, we apply the X-ray virtual telescope PHOX 
to generate synthetic observations of the simulated clusters with both Chandra and 
ATHENA, the latter used as an example for the performance of very high-resolution 
X-ray telescopes. From Chandra spectra we extract global properties, e.g. luminosity 
and temperature, and we accurately estimate the gas velocity dispersion along the 
line of sight achievable from the broadening of emission lines from heavy ions (e.g. Fe) 
resolved in ATHENA spectra. Given the good agreement found between simulations 
(true, intrinsic solution) and mock observations (detectable amplitude of non-thermal 
velocities), we further extend the analysis to the relation between non-thermal ve- 
locity dispersion of the gas and the Lx — T scaling law for the simulated clusters. 
Interestingly, we find a clear dependence of slope and scatter on the selection criterion 
for the clusters, based on the level of significance of non-thermal motions. Namely, 
the scatter in the relation is significantly reduced by the exclusion of the clusters, for 
which we estimate the highest turbulent velocities. Such velocity diagnostics appears 
therefore as a promising independent way to identify disturbed clusters, in addition 
to the commonly used morphological inspection. 
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1 INTRODUCTION 

During the mass assembly in galaxy clusters, interactions 
between sub-haloes and merging events can generate sub- 
stantial streaming motions and turbulence in the hot gas 
filling the cluster potential well. Additionally, a number of 
processes taking place in galaxy clusters, especially in the 
inner regions, are most likely responsible for the transfer of 
energy from large modes into smaller modes, causing rota- 
tion, streaming and, mainly, turbulent motions to establish 
in the intra-cluster medium (ICM). Among these physical 
processes, mergers and sloshing of dark matter cuspy cores 
are believed to cause large scale motions, and, in addition, 
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the AGN activity and its interaction with the surrounding 
gas can cause turbulence in the central region. 

Numerically, the non-thermal components of ICM mo- 
tions have been investigated by means of hydrodynam- 
ical simulations of galaxy clusters, which uniquely pro- 
vide full 3D information on the gas velocity field. The es- 
tablishment of ICM bulk, streaming and rotational mo- 
tions during the growth and assembly of simulated galaxy 
clusters is believed to contribute, especially in the cen- 
tral part of the system, to the cluster pressure support 
and therefore to the virial estimate of the total mas s up 
to significant fractions (e.g. Pawl et al. 2005; Fang et al.l 
l2009l : lLau et al.ll2009l . l201ll : lBiffi et al.ll2011h . Moreover, sev- 
eral studies on SPH and AMR simulations of cluster-like 
haloes have addressed the tu rbulent velocity field in clus- 
ters (e.g., see recent work bv lPaul et al.|[201ll : IVazza et al.l 
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I2OIII : IValdarninil l201lh . estimating the pressure support 
due to chaotic motion s to be of order of ^ 20 — 30% of 
the t otal pressure (e.g. Norman & Bryan 1999; Dolag et al. 
20051: fla pichino & Niemeyer 2008; lapichino et al. 2011; 



Vazza et al.ll2009lJ201ll V Tighter constraints on the gas ve- 



locity field are necessary in order to obtain precise mea- 
surements of the total gravitating mass, which is the most 
important, intrinsic quantity to determine. In particular, the 
account for non-thermal motions is essential for mass esti- 
mates based on X-ray global properties (e.g. gas density and 
tem perature), via the assumption of hydrostatic equi librium 
(e.g. iRasia et al .11 20061 . 12OI2I : IPiffaretti k Valdarnini| [2008) . 

The presence of non-thermal motions in the gas ve- 
locity field within galaxy clusters is also suggested by sev- 
eral observational evidences, coming from radio observations 
of p olarized synchrotron emission in cluster radio g alaxies 
(e.g. ICassano k, Brunettill2005l : iBonafede et al.l 2010lV mea- 
surements of the resona nt scattering effects (jChurazov et all 
I2OIOI : [Zhurayleva et al. 201Q) and study of the fluctuations 
i n pressure (ISchuecker et^tll |2004| ) and surface brightness 
(|Ch urazov e t al.l l201lh maps, obtained with X-ray tele- 
scopes. However, mainly indirect indications of ICM turbu- 
lence have been possible so far. Only with the XMM-Newton 
RGS spectrometer weak upper limits have been set on the 
turbulent v elocities in a set of g alaxy clusters, as recently 
discussed in I Sanders et al.1 (l201lh . 

Future X-ray missions, like ASTRO-H, will allow us to 
achieve direct estimations of the ICM non-thermal velocities 
with great accuracy, thanks to the high spectral resolution 
expected to be reached. In fact. X-ray high-precision spec- 
troscopy potentially offers one of the most promising ways 
to directly measure such gas motions, detectable from the 
detailed study of the shape and centroid of resolved spectral 
emission lines. 

Theoretically, the non-thermal component of the gas 
velocity, along the line of sight (l.o.s.), can be very well 
constrained by studying the shape of heavy-ion emission 
lines in the X-ray spectra, for which the broadening can 
be significantly more sensitive to non-thermal veloc ities of 
the g as rather than to thermal motion (e.g. Rebusc o et al.l 
I2OO8I ) . The expectations for such line diagnostics are related 
in particular to the most prominent emission line in X-ray 
spectra, namely the ^6.7 keV line from helium-like iron. In 
fact, the large atomic mass of the FeXXV ion significantly 
reduces the thermal line broadening and the line width turns 
out to be definitely more sensitive to turbulent g as motions 
(inogamov Sunvaevlliooi : ISunvaev etal l liooi ). 

Here we discuss the prospect of using high resolution 
spectra to detect the amplitude of non-thermal gas motions 
with the aim of characterising more precisely galaxy clusters 
and observed relations between their X-ray properties, like 
the Lx — T scaling relation. 

To this scope, we employ a set of numerically simulated 
clusters, obtained with the TreePM/SPH parallel code P- 
G AD GETS including a large variety of physical processes 
to describe the baryonic component, and perf orm X-ray 
synth etic observations of the haloes with PHOX (|Biffi et al.l 
I2OI2I ). The paper is structured as follows: first, we will de- 
scribe the simulated dataset of galaxy clusters (Section [2]) 
and the generation of mock X-ray spectra with the Chan- 
dra telescope and the X-ray spectrometer originally planned 
for the ATHENA satellite, as prototype for X-ray spec- 



troscopy at very high energy resolution. In Section [S] we 
will describe the analysis performed to obtain global prop- 
erties, such as luminosity and temperature, from the Chan- 
dra spectra and to estimate the gas non-thermal velocity 
dispersion, from the velocity broadening of the iron lines in 
the high-resolution ATHENA spectra. The analysis on the 
gas velocity dispersion, calculated directly from the simu- 
lation, is presented in Section 14.11 The comparison against 
the detectable velocities obtained from the synthetic high- 
resolution spectra is then discussed in Section 221 Using the 
Chandra mock observations, we explore the Lx — T scaling 
relation for the simulated clusters in Section 14.31 In Sec- 
tion [5] we discuss the effects on the best-fit relation, due to 
the non-thermal fraction of the ICM motions, and the rela- 
tion between velocity diagnostics and cluster internal state. 
Our conclusions are finally summarized in Section (6] 



2 THE SAMPLE OF SIMULATED CLUSTERS 

The sample of simulated cluster-like haloes analysed here 
has been extracted from a cosmological, hydrodynami- 
cal simulation performed with the TreePM/SPH, paral- 
lel code P-GADGET3. In this extended ve rsion of GAD- 
GET2 tepringel et al. 1 12OOII : 'Springellliooi), a vast range 
of baryonic physics, at a high level of detail, is included, 
such as cooling, star fo rmation and supernova-driven winds 
(jSpringel HernquistI 12003 ) , chemical e nrichment from 
stellar pop ulation, AGB stars and SNe (iTornatore et aP 
I2OO4I . l2007l \ low-viscosity SPH (Dolag et al.""2005'), black- 
hole growth a nd feedback from AGN (Springel et al. 20051: 

iFabian et al. 1 [mfl). 

The simulated box has a side of 
352/i~^ Mpc, in comoving units, resolved with 2x594^ parti- 
cles, which results in a mass resolution of ttidm = 1.3 x 10^° 
and rrigas — 5.2 x 10^, for dark matter (DM) and gas 
particles, respectively. For the simulation, and throughout 
the following, the cosmology assumed re fers to the 7-year 
WMAP estimates (jKomatsu et al.ll201lh . i.e. = 0.268, 
= 0.728, = 0.044, erg = 0.776 and h = 0.704. 
The sample of clusters consists of 43 objects, selected 
among the most massive haloes in the simulated box, for 
a snapshot of the simulation at z — 0.213. The selection 
criterion adopted requires the cluster mass encompassed by 
i?50cEl, M500, to be > 3 X 10^^ M©, at the redshift con- 
sidered. 



3 X-RAY SYNTHETIC OBSERVATIONS 

The mock X-ray observations have been performed for the 
selected haloes of the sample by means of the virtual X-ray 
simulator PHOX. For a detailed d escription o f the m ethod 
implemented in PHOX, we refer to lBiffi et aP (|2012l ). 



^ Note that i^soo is defined here as the radius enclosing the region 
of the cluster whose mean density is 500 times the mean density 
of the Universe. This encompasses therefore a larger region with 
respect to the usual definition of i^5oo, where the overdensity is 
instead defined with respect to the critical density of the Universe. 
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3.1 Generation of the virtual photon lists 

The snapshot of the hydro dynamical simulation considered, 
referring to redshift z = 0.213, has been first processed with 
Unit 1 of PHOX as a whole, generating the virtual photon 
cube associated to the X-ray emitting gas component in the 
simulated cosmic volume. 

For all the gas particles, model spectra have been cal- 
culated and sampled with packages of photons, each of 
them in the restframe of the corresponding emitting par- 
ticle. The model adopted to calculate each theoretical spec- 
trum was an abso rbed, single-temperature VAPEcEl mode l 
(jSmith et al.ll200lh . implemented in XSPEcQ (|Arnaudlll996h . 
Here, we make use of the parameter within XSPEC (v. 12) 
to include the thermal broadening of emission lines in the 
model spectra. Temperature, density and chemical abun- 
dances (see Appendix [Al for details on the implementation 
of metal composition in PHOX) are obtained directly from 
the hydro-simulation output. Additionally, we fixed the red- 
shift at the value of the simulated data cube (z — 0.213) 
and the equivalen t hydrogen column density for the WABS 
absorption model ([Morrison McCammonll 19831 ) to a fidu- 
cial value of 10^° cm~^. At the end of this first stage, the 
photon database associated to the simulation output con- 
tains roughly 10^ photon packages 10^ photons, in to- 
tal), for fiducial values of collecting area and exposure time 
of 1000 cm^ and IMs, respectively. 

With PHOX Unit 2 we then assume to observe the pho- 
tons from a line of sight (l.o.s.) aligned with the z-axis of the 
simulation cube and correct photon energies for the Doppler 
Shift due to the emitting-particle motions along this l.o.s.. 
During this geometrical stage of the process, we also select 
cylindrical sub-regions along the l.o.s., corresponding to the 
43 selected cluster-like haloes. For each cluster halo, the 
selected sub-region is centered on the centre of mass and 
encloses the region within i^soo, in the xy plan, throughout 
the box depth (along the 2;-axis). For the time frame consid- 
ered here and the cosmology used, the luminosity distance 
between the observer, positioned along the positive z-axis, 
and the observed region is 1047.6 Mpc, in physical units. 
Typically, we obtain 1 — 2 x 10^ photons per halo, for lOOks 
exposure. 

The photon lists produced in this general way, are con- 
volved by PHOX with real instrumental responses of Chan- 
dra and also for a high-energy resolution response, origi- 
nally planned for the X-ray spectrometer of ATHENA. The 
synthetic spect ra are fitted b y means of the X-ray package 
XSPEC v.12.6.0 (|Arnaudlll99"6l ). 

3.2 Chandra synthetic spectra 

To create Chandra synthetic spectra, we use the ARF and 
RMF of the Chandra ACIS-I3 detector aimpoint. The FoV 
of Chandra, which is 17' x 17^ corresponds to a physical 
scale of 3.52 Mpc per side, for the given cosmology and red- 
shift. This encloses typically the region within 1 — 1.2i?5oo 
for most of the clusters in the sample, except for the 7 most- 
massive haloes, for which the Chandra FoV captures a re- 
gion slightly smaller than the one out to i?5oo- In such cases. 



See |http: / /heasarc.gsfc. nasa.gov/xanadu/xspec/manual/XSmodelApec. htm l 
See |http: / /heasarc.gsfc.nasa.gov/xanadu/xspec/| 



however, one can in principle compose a mosaic with mul- 
tiple Chandra pointings in order to cover the clusters up to 
i?500- Therefore, we convolve the photon lists corresponding 
to the whole i?5oo region with the response of Chandra and 
analyse the corresponding spectra, for all the clusters in the 
sample. 

3.2.1 Temperature and bolometric luminosity 

We fit Chandra spectra either with an absorbed, single- 
temperature or with a two-temperature APEC model, de- 
pending on the goodness of the single-temperature fit. The 
spectra are re-grouped requiring a minimum of 10 counts per 
energy bin. Whenever the spectral fit was still poor, spectra 
were instead re-grouped with a minimum of 20 counts. The 
absorption is fixed to the fiducial value adopted to generate 
photons. Given the characteristics of the Chandra response, 
we assumed the redshift to be fixed at the initial value of 
the simulation and only temperature, metallicity and nor- 
malization were free in the fit. 

From the spectral best fit in the 0.5 — 10 keV range we 
directly obtain the temperature of the clustefl 

The total X-ray luminosity, Lx, is extrapolated from 
the best-fit spectral model out to the entire energy range 
defined by the ACIS-I3 detector response, to obtain an ap- 
proximate bolometric luminosity. 

3.3 High— resolution synthetic spectra 

In the present work, we consider the ideal instrumental 
response of the spectrometer designed for the ATHENA 
satellite, which has, however, not been chosen as the next 
ESA L-class mission in the recent mission selection pro- 
cess. Nevertheless, we use this as proptotype and goal for 
next-generation high-resolution X-ray spectrometers (sim- 
ilar to ASTRO-H), fundamental to explore the details of 
hot plasma emission. 

Therefore, the mock spectra are simulated from the 
ideal photon lists of the 43 haloes in the sample by using 
the latest RSP response matrix planned for the ATHENA 
X-ray Microcalorimeter Spectrometer (XMS). For the given 
cosmological parameters and redshift, the small FoV of such 
instrument (2.4' x 2.4') encloses a region of 497.52 kpc per 
side. We center the field of view on the centre of the cluster, 
probing therefore the very central part. 

3.3.1 Velocity broadening 

Given the high energy resolution provided by the XMS 
spectrometer, spectra emission lines are very well resolved. 
In fact, we restrict the spectral analysis to the observed 
5 — 6.5 keV energy band, containing the iron Ka complex 
around 6.7 keV (rest-frame energy), and fit with an ab- 
sorbed BAPEc[fl model (see Fig.[T] for an example with one 
of the presented haloes), i.e. a velocity- and thermally- 
broadened emission spectrum for collisionally-ionised diffuse 

^ In the case of haloes fitted with a two-component model, we 
assume here the hotter temperature to be representative of the 
dominant gas component. 

^ See http: / /heasarc. nasa.gov/xanadu/xspec/manual/XSmodelBapec. html] 
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Figure 1. Zoom onto the 5 — 6.5 keV energy band, containing 
the He- and H-hke iron emission hnes (at the rest-frame energy 
- 6.7 keV and - 6.96 keV), from the synthetic ATHENA XMS 
spectrum of one of the sample cluster. In the figure we show the 
best-fit BAPEC model to the data (green curve). 

gas. The model assumes the distribution of the gas non- 
thermal velocity along the l.o.s. to be Gaussian and the ve- 
locity broadening is quantified by the standard deviation, a, 
of this distribution. 



4 RESULTS 

4.1 Velocity statistics 

Our first goal is to investigate the properties of the ICM 
velocity field in the simulated data directly, in order to es- 
tablish the level of reliability of the mock data results, and 
therefore the validity of the expectations for high-resolution 
spectroscopy. Gas particles in the simulated haloes have 
been selected to reside within i?5oo in the xy plane, exclud- 
ing the low-temperature {Tgas < 3 x 10^ K) and star form- 
ing (gas particles with > 10% of cold fraction) phase of the 
gas. This selection is done to resemble in the most faithful 
way possible the X-ray emitting gas residing in the region 
observed with the X-ray virtual telescope PHOX. In the 
following, we will always refer to regions in the xy plane, 
considering the cylinder-like volume along the line of sight, 
aligned with the ^-axis of the simulated box. 

By definition, we calculate the weighted value of the 
standard deviation of the gas velocities along the 1 .O.S., (Tuo 
as 

= — — ' 

where w represents the quantity used to weight the veloc- 
ity (e.g. the mass or the emission measure, EM, of the gas 
particle), Vi is the l.o.s. component of the velocity for each 
particle and Vw is the weighted mean value of the l.o.s. ve- 
locity. 

From Fig.O the relation between the value of cr5oo,m 
(i.e. weighted by the particle mass) and M500 shows that, 
in general, more significant motions of the gas are expected 
in more massive systems. Ideally, the mass-weighted values 
should trace in the most faithful way the intrinsic quanti- 
ties characterising the cluster, since closely related to the 



potential well of the halo. However, X-ray observations of 
the ICM are more sensitive to the emission properties of 
the gas, like the emission measure (EMjl, and rather pro- 
vide estimates for an EM- weighted-like velocity dispersion. 
We therefore investigate the relation between mass-weighted 
and EM-weighted velocity dispersion, (7500, m and (t^qq^em 
respectively. 

The comparison is shown in Fig.[3l wherein both values 
of (7500, are calculated for the gas residing within i?5oo, 
in the plane perpendicular to the l.o.s.. The relation found 
between the two definitions of (t^qq^w is not coincident with 
the one-to-one relation (overplotted in red in the figure), 
as the EM-weighted value is likely to be affected by the 
thermo-dynamical status of the gas, that is by processes 
such as turbulence, merging and substructures. Despite this, 
the two values are fairly well correlated. 

For the purpose of our following analysis and the com- 
parison against synthetic X-ray data, however, we decide 
to use the EM-weighted velocity dispersion, (t^qq^em^ which 
is more directly related to the X-ray emission of the gas, 
because of the proportionality between the normalization of 
the X-ray spectrum and the gas EM itself. 

In order to probe the global^ dynamical structure of the 
ICM we would need to observationally measure the gas ve- 
locity dispersion within the whole i?5oo region. Therefore, 
we explore the relation between the estimated value of the 
velocity broadening along the line of sight in different re- 
gions of the cluster, shown in Fig.^J It is evident from the 
figure that the value calculated for the gas within i?5oo cor- 
relates linearly with the value computed in smaller, internal 
regions, namely for r < 0.3i?5oo (upper panel) and for the 
region covered by the FoV of ATHENA, - 0.15i?5oo (lower 
panel) . With respect to the one-to-one correlation (red line 
in Fig. 31), however, outliers are present in this sample, show- 
ing that prominent substructures in the velocity field of the 
gas must be present in the observed regions around the clus- 
ters. Here, as in Fig. [21 deviations can also be related to sub- 
haloes residing in the projected i?500, i-e. in the region ob- 
served along the line of sight, but not necessarily comprised 
within the three-dimensional i^soo- 

Therefore, the level of complexity in the spatial struc- 
ture of the ICM velocity field can be singled out by the 
comparison between the velocity dispersion calculated in 
the i?500 region and the values corresponding to smaller, 
inner regions. 

Nevertheless, the relations discussed ensure that, for the pur- 
poses of this work, we can safely: 

(i) assume the EM-weighted velocity dispersion instead of 
the mass-weighted value to trace the intrinsic velocity struc- 
ture; 

(ii) focus on the expected value for the whole R^m region of 
the cluster. This traces fairly well the internal motions also 
within smaller regions, as those covered by instruments with 
a smaller FoV, with the exception of those cases where self- 
bound substructures in the outskirts, or more likely along 
the line of sight, are clearly present and might be important 
to consider. 



We recall here that EM = / neUndV . 
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Figure 2. Theoretical value of the mass-weighted velocity disper- 
sion, am, calculated within R500 (in the plane perpendicular to 
the 1.0. s. direction), reported as function of the halo mass M500, 
in Mq. 
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Figure 3. Theoretical value of the EM-weighted velocity disper- 
sion, (7500, £;m versus the mass-weighted value, (7500, m, in km/s. 
Both values are calculated for the region within i^5oo 5 in the plane 
perpendicular to the l.o.s. direction. Overplotted in red is the 
curve referring to the one-to-one relation. 



4.2 Comparison against synthetic data 

The ICM velocity dispersion calculated directly from the 
simulation can here be used to compare against the mock 
ATHENA data, which is the prototype for the high- 
resolution spectroscopy (X-ray microcalorimeter spectrom- 
eter, XMS) required to measure the gas velocity dispersion 
along the line of sight from resolved emission lines. 

This will allow us to (i) test whether predictions from 
simulations are indeed comparable to observational data 
and, in turn, (ii) constrain the reliability of high-resolution 
X-ray spectroscopy to derive direct measurements of the 
ICM velocity field. 

Fig. [5] shows the comparison between expectations pro- 
vided by the simulated data (black diamonds) and results 
from analysis of the synthetic ATHENA spectra (blue aster- 
isks) . 
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Figure 4. Relation between the EM-weighted velocity disper- 
sion within cr5oo,£'M: cind the analogous values calculated 
for: (top panel) the region within 0.3i^5oo and (bottom panel) the 
region covered by the FoV of ATHENA, (i.e. ~ 0.15i?50o). Over- 
plotted in red is the curve referring to the one-to-one relation. 



The expected velocity dispersion of the gas particles re- 
siding in the region corresponding to the ATHENA FoV is 
calculated according to Eq.[Tl weighted by the EM, while 
the values derived from the X-ray spectra are obtained as 
described in Section [3.3.1l with error bars corresponding to 
the la errors to the best-fit values. As shown in the Figure, 
we find very good agreement between simulation (intrinsic, 
"true" solution) and synthetic spectral data (observational 
detections), namely for ^ 74% of the haloes the spectral 
analysis of the iron lines provides a measure of the gas ve- 
locity dispersion, along the l.o.s., within 20% from the ex- 
pected value (purple, shaded area). We find, in particular, 
that ^ 50% of the halos show agreement at a level bet- 
ter than ^ 10% (internal, pink shaded area in Fig.[5|). We 
remark here that the reference number, or halo id in the Fig- 
ure, is associated to the haloes of the sample in an increasing 
order for increasing Msoo- 

The deviation between expected and measured velocity 
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Figure 5. Comparison between the theoretical expectation of 
the velocity dispersion (Jv^ calculated directly from the simulation 
(black diamonds and shaded areas), and the value obtained from 
the spectral fitting of the synthetic XMS spectra obtained with 
PHOX (blue asterisks with error bars). The id numbers of the 
43 haloes in the sample (x-axis) are ordered according to the 
increasing halo mass, Msoo- 




Figure 6. Deviation of the best-fit c^^, from the expected value 
for all the 43 haloes in the sample. 



dispersion, referring to the XMS FoV as in Fig.[5l has been 
quantified as 



(2) 

and its distribution for the sample is reported in Fig.[6l 
Clearly, the distribution of 5 is peaked around the zero value, 
reflecting the very good agreement between the two mea- 
surements. 

However, we also find outliers in the sample that show 
deviations up to ^ 43%. 

4.2.1 Extreme cases in the sample 

Given the distribution of the deviations (Fig. [6]), we focus 
onto two subsets of haloes in the sample for which the devi- 
ation between simulation and mock data is very minor and 
most prominent, respectively. 

In Fig. [71 we show the l.o.s. velocity structure as cal- 
culated from the simulation directly, for both subsets. In 



particular, the 5 haloes for which mock and simulated data 
most disagree are shown in the first and second columns, 
while the least-deviating haloes are presented on the right- 
hand-side of the Figure (third and fourth columns). The 
rows in the Figure correspond to a decreasing level of devi- 
ation, or agreement, from top to bottom. 

Histograms in the central columns of the Figure, report 
the EM distribution as function of the l.o.s. velocity, for the 
gas particles in the ATHENA FoV. Black curves refer to all 
the gas in the region, while the overplotted red histograms 
only account for the hot-phase gas, i.e. particles for which 
kT > 2 keV. The reason for selecting the hot gas is that it 
mostly contributes to the iron line emission, from which the 
velocity broadening is measured. 

It is clear that in the most-deviating cases there are sub- 
stantial substructures within the gas velocity field. The red 
histograms for the haloes that show best agreement (right- 
hand-side column), instead, reflect more regular distribu- 
tions of the EM as function of the vi.o.s., indicating more 
regular velocity fields. 

Observationally, the value estimated from the broad- 
ening of the spectral lines is assumed to be the dispersion 
of the Gaussian distribution that best fits the line shape. 
Therefore, a more detailed comparison should involve the 
dispersion of the Gaussian function matching the (red) dis- 
tribution shown, instead of the theoretical value calculated 
as in Eq.[T] The green curves in Fig. [7] define the Gaussian 
fits to the red distributions, whose cr^^^** is more directly 
comparable to the mock spectral results. 

As an additional comparison, we also overplot the Gaus- 
sian curve (blue asterisks) constructed from the theoreti- 
cal estimation of the EM-weighted values for gas velocity 
dispersion (Eq.[T]) and mean l.o.s. velocity. In the most- 
deviating clusters, the evident differences among the green 
and blue Gaussian curves effectively reflects the deviations 
discussed above (see, e.g., the halo with the largest devia- 
tion, top row, left panels in Fig.[71). In particular, the green 
best-flt Gaussian clearly fails to capture all the features of 
the multi-component velocity distribution, as it most likely 
happens during the spectral flt. The low level of deviation 
found for the "best" halo set (right-hand-side columns) is 
indeed shown by the good match between best-fit (green 
curve) and theoretical (blue asterisks) Gaussian overplotted 
to the EjM-vi,o.s. distributions. 

A visualization of the thermo-dynamical structure of 
these clusters is also shown in the first and fourth columns 
of Fig. [71 for the subsets of most-deviating and best haloes, 
respectively. The maps show the vi,o.s. — kT distribution for 
the gas residing within the selected FoV, color-coded by 
EM. 

Especially for the most deviating haloes, the sub- 
structures in the velocity field unveiled by the histograms 
are clearly visible, combined with the gas thermal structure. 

Despite the deviations discussed, a good overall agree- 
ment is found between the intrinsic amplitude of the gas 
velocity dispersion and the velocity broadening measured 
directly from mock spectra. This confirms the promising po- 
tential of such well-resolved observations, which would most 
likely allow us to derive reliable and precise constraints on 
the ICM velocity field. 

Moreover, given that the simulation analysis suggests 
the gas velocity structure in the innermost region (e.g. that 
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5 MOST-DEVIATING HALOES 5 LEAST-DEVIATING HALOES 




Figure 7. Comparison among the five most-deviating (left) and least-deviating (right) haloes, according to Fig. [6] for the central region 
corresponding to the 2A' x 2A' FoV. The level of deviation, or agreement, decreases from top to bottom. From left to right, the columns 
refer to: (1) 'u^o.s. ~ distribution for the gas residing within the selected FoV, color-coded by EM, for the most-deviating haloes; 
(2) distribution of EM as function of gas l.o.s. velocity, for the five most-deviating haloes; (3) same as (2) and (4) same as (1), but 
considering the least-deviating haloes. In the histograms shown in columns (2) and (3), the curves refer to: all the gas particles (black); 
hot-phase gas {kT > 2 keV, red); Gaussian best-fit to the hot-gas distribution (green); theoretical Gaussian distribution reconstructed 
from the estimated dispersion of the hot-gas distribution (blue asterisks); 



covered by the FoV of ATHENA) to be closely traced by 
that within i?5oo (see Section l4.ip , we will use the latter 
for our further investigation. Eventually, limitations due to 
small fields of view could be observationally overcome by 
multiple pointings covering a larger region. 



4.3 Lx — T scaling relation 

Here, we investigate the impact of the ICM velocity struc- 
ture on X-ray global properties by focusing on the Lx — T 
relation. The main motivation behind this choice is that, 
on one hand, the luminosity Lx, is very well measured in 
X-ray surveys (e.g. with Chandra, XMM-Newton, or the 
up-coming eRostita), and, on the other hand, the tempera- 
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Figure 8. L — T scaling relation from Chandra mock spectra 
(green asterisks), for the region of the clusters encompassing i^soo- 
Overplotted are the results from P09 (black diamonds). Best-fit 
relations are also overplotted (solid lines for P09; dotted lines for 
simulated haloes of the present work) and correspond to the two 
linear regression fitting methods adopted: BCES (L|T) (red) and 
BCES Orthogonal (blue). 



ture T provides a good mass proxy, since it is tightly related 
to the total gravitating mass, which is the most fundamental 
quantity to characterise a cluster. 

Lx, also denoted as "bolometric X-ray luminosity", is 
usually the X-ray luminosity extrapolated to the whole X- 
ray band, 0.1 — 100 keV, instead of being calculated in a 
narrow energy band. Nonetheless, in our case, the computa- 
tion of the luminosity is limited to the largest energy band 
defined by the ACIS-I3 response matrix (i.e. 0.26 — 12 keV), 
since the difference introduced with respect to the expected 
bolometric X-ray luminosity is found to be minor. 

In Fig. [8] we show the Lx — T relation calculated for 
the simulated haloes (asterisks), from the Chandra synthetic 
observations of the R500 regions through the l.o.s., obtained 
with PHO X. As a comparis on, we also show the data pre- 
sented in (|Pratt et al.1 l2009l . also, P09) for a sample of 31 
nearby galaxy clusters (black diamonds), selected only in 
X-ray luminosity from the Representative XMM-Newton 
Cluster Structure Survey (REXCESS). 



4-3.1 Fitting procedure 

The Lx — T relation was fitted with a power-law functional 
form, that is 



h-'(z)Lx=C(T/nT, 



(3) 



where Lx is given in units of 10 erg/s. To = 5 keV and 
the assumed cosmology is accounted for in h^{z) — VLq{1 + 
zY -\- Qa. The fit was performed using linear regression in 

the log-log plaim 

As done in Pratt et al. (2009"), as well as in several sim- 
ilar works (e.g. Reiprich & Bohringer 2002; Arnauc Tet al.l 
I2OO5I : iMittal et al.ll201lh . we adopted the BCES (Bivariate 
Co rrelated Errors and i ntrins ic Scatter) regression method 
bv lAkritas Bershadvl (|l996l ). which accounts for both the 
errors in Lx and in T, as well as the intrinsic scatter in the 



data. The BCES algorithm allowed us to calculate the best- 
fit values of slope [a) and normalization (C) for four differ- 
ent regression methods, amongst which we restrict our atten- 
tion to the BCES (L|T) and the BCES Orthogonal methods. 
Our primary goal is to find the best fit which minimizes the 
residuals of both variables at the same time, orthogonally to 
the linear relation. This is given by the BCES Orthogonal 
method. Additionally, we also explore the resu lts given by 
the BCES (L|T) fitting method (analogously to lPratt et ahl 
(j2009)), which minimises the residuals in Lx- Reasons for 
this rely on the fact that Lx can be treated as the dependent 
variable, while the temperature can be assumed to be the 
"independent" one, as closely related to the cluster mass, 
which is the intrinsic quantity characterizing the system. 
Given the statistical uncertainties on both variables, 
and (JXi, the raw scatter was estimated using the or- 
thogonal distances to the regression line, weighted by the 
errors. Namely, for a linear relation of the form Yi — aXi-\-P 
in the log-log space, the raw scatter is 



2 _ N/{N - 2) 



(4) 



where N is the number of data points in the sample, Yi 
log(Lx,0, y^ = log(TO, P = logC, and 



2,22 



(5) 



The intrinsic scatter was estimated from the difference be- 
tween the raw and the statistical scatter, in quadrature 

/ 2 _ 2 I 2 \ 

y^raw ^ stat intrinsic ) • 

Best-fit relations are overplotted in Fig. [8] for both 
observational data (solid red and blue lines, for the 
BCES (LI T) and the BCES Orthogonal method respectively, 
from Pr att et al.|[2009 ) and for Chandra synthetic observa- 
tions of the simulated sample (dotted lines). The linear re- 
lations found for the simulated clusters are overall shallower 
than the observed ones, and, among the two fitting methods 
considered, the BCES (L|T) method still provides a shal- 
lower slope than the BCES Orthogonal case. 

In particular, we find a slope a(L|T) = 1-97 =b 0.23 
for the BCES (L|T) fit and aortho = 2.78 ± 0.3 for the 
BCES Orthogonal fit. As for the normalization of the best- 
fit relations, we find C^L\T) = 6.81 ± 0.39(10^^ erg/s) and 
Cortho = 6.11 ± 0.34(10^^ erg/s), respectively. 

We note, however, that our cluster sample probes a 
smaller dynamical range with respect to observations and, 
in particular, lacks low-temperature haloes, whose presence 
might contribute providing tighter constraints on the slope 
of the relation. 



5 DISCUSSION 

5.1 Lx — T relation: effects of the velocity 
structure 

In Fig.[9]we include an additional parameter to further char- 
acterize the cluster state. In particular we show the Lx — T 

^ In the log space, errors are transformed as Alogo: = loge x 
(Ax)/ (2a::), where Ax is the difference between the upper and 
lower boundary of the error range around the quantity. 
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Figure 9. Lx — T scaling relation from Chandra mock spectra, for the region of the clusters encompassing R500. Color code: value of 
fj, = av(< R500) / thermal ^ Calculated from gas velocities in the region encompassed by R500. Observational data points (black diamonds) 
are from P09, with cool core clusters marked with thicker black diamonds. 



relation for the simulated haloes, as in Fig. [51 color-coding 
the data points according to the increasing ratio of 



/i = crv{< R50o)/crth, 



ermal 5 



(6) 



calculated for the region encompassed by i?500- 

The velocity dispersion ay corresponds to the EM- 
weighted value calculated directly from the gas particles in 
the simulation (see Section [4.1|) . The thermal velocity dis- 
persion, athermah IS the cxpcctcd three-dimeusioual value 
for the ICM temperature T, reported in the x-axis of the 
relation. To this purpose of normalizing the non-thermal 
velocity dispersion to a value characteristic for the halo po- 
tential well, the choice of the three-dimensional thermal ve- 
locity is equivalent to the ID value, apart from the constant 
scale factor v^- Small values of fi indicate a low level of 
non-thermal velocity with respect to the characteristic ther- 
mal velocity dispersion of the gafl For the extreme case of 
/i ^ 1, the non-thermal velocity dispersion would equal the 
thermal value. 

The distribution of /i for the 43 simulated clusters is 
shown in Fig. 1101 A significant number of haloes 35%) 



^ In fact, the thermal component of the velocity is not included 
within ay , which traces macroscopic motions of the gas elements 
in the simulation. Analogously, the value measured from the ve- 
locity broadening of spectral lines, used for comparison in Sec- 
tion ^2] does not include the thermal component either. 
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Figure 10. Distribution of /i = (Jv(< R500) / o- therm ah for the 
sample of 43 simulated clusters. Green, dashed lines correspond to 
the values /J^max = 0.5, 0.35 used to identify the two sub-samples. 



are characterised by a gas velocity dispersion which is larger 

than OAathermal- 

In particular, by examining the distribution of fi within 
the sample, we decide to extract two additional sub-samples 
from the 43 haloes, selected to have a maximum of /Umax = 
0.5, 0.35, respectively. The first sub-sample is intended to 
exclude the most prominent outliers in the /x distribution. 
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shown in Fig. 1111 where we report both the BCES (L|T) (red 
squares) and the BCES Orthogonal (blue squares) values. 

For all the quantities studied, the best-fit values 
generally increase with the increasing amplitude of velocity 
dispersion, relative to the characteristic thermal value. 
This has significant implications for the intrinsic scatter 
(calculated here for the luminosity, (T\rvL,intrinsi^- In fact, 
the trend found reflects how the scatter in luminosity of the 
Lx — T scaling relation is sensitive to the baryonic physics 
and can be closely related to complex or disturbed config- 
urations of the ICM velocity field, which are quantified by 
large values of /x. Essentially, we find that the introduction 
of clusters with significant non-thermal velocity dispersion, 
with respect to their typical thermal velocity, augments the 
scatter in the sample about the best-fit Lx — T relation. 
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Figure 11. Dependence of normalization (C, top), slope (o;, mid- 
dle) and intrinsic scatter {(J\nL, intrinsic-, bottom) of the L — T 
relation on the maximum /i used to select the corresponding sub- 
sample of clusters. 



for which the velocity dispersion ay exceeds 50% of cr thermal- 
The smaller sub-sample contains the haloes with the small- 
est fraction of /i, indicating that their expected thermal ve- 
locity is dominant with respect to cr^. 

In order to quantify the effects of l.o.s. velocity struc- 
ture on the resulting scaling relation, we investigated the 
dependences of best-fit slope, normalization and scatter on 
the value of /i, by fitting the Lx — T relation for the original 
sample and for the two sub-samples, separately. Results are 



We note from the relation shown in Fig. (9] that the ad- 
ditional characterization of the cluster through the line-of- 
sight non-thermal velocity dispersion defines a different pic- 
ture with respect to observations. The observed behaviour 
of cool-core and disturbed galaxy clusters shows in fact a 
distinct separation of the two populations in the Lx—T rela- 
tion (s ee data points reported in the Figure from lPratt et al.l 
(|2009h V where cool-core clusters are generally found to oc- 
cupy the upper envelope of the best-fit relation (black, thick 
diamonds in Fig. [9]), while morphologically disturbed clus- 
ters mostly reside in the lower one. In Fig.[9l however, we find 
that haloes with the most significant degree of non-thermal 
motions populate mostly the region above the best-fit curve. 
According to this velocity diagnostics, these haloes might be 
classified as disturbed. We interpret this apparent inconsis- 
tency with observations as a different probe for deviation 
from the regularity of the haloes. Indeed, a deeper analysis 
of the simulated sample indicates that all the clusters have 
central electron densities which are sufficiently low to be 
identified as non-cool core clusters, which remove the dis- 
crepancy with observed clusters entirely, as further discussed 
in Section 15.31 Namely, within the population represented 
by our sample, the quantification of non-thermal velocity 
dispersion of the gas along the line of sight consitutes an 
additional, complementary method to further discriminate 
the halo dynamical state. 



5.2 Velocity diagnostics: prominent outliers 

The velocity diagnostics performed and the distribution of 
the 11 value clearly single out some haloes in the high- 
velocity tail. This subsample of 5 clusters includes the most 
prominent outliers in the Lx — T relation, and the exclusion 
of these haloes turns out to reduce the scatter in the best-fit 
relation. Therefore, we further discuss here their thermo- 
dynamical structure from the simulation directly, in order 
to provide an interpretation of the high gas velocities found. 

In Fig. [12] we show the spatial distribution of tempera- 
ture (kT), l.o.s. component of the gas velocity {vi,o.s.) and 
emission measure (EM) of the 5 outliers, as well as the 



^ We report the final value by considerin g the natural logarithm, 
In, for comparisons to lPratt et al 1 (l2009h . although the equations 
and the logarithm space mentioned in the paper usually refer to 



log = logi 
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Figure 12. Maps of the ICM thermo-dynamical status for the five most-significant out hers singled out by the velocity diagnostics. 
Columns, from left to right: (1) vi^o.s. — kT map, color-coded by EM (considering gas particles within i^soo); (2) EM-weighted, tem- 
perature map, projected along the l.o.s.; (3) EM-weighted, l.o.s. -velocity map projected onto the xy plan (i.e. along the l.o.s.); (4) EM 
map projected along the l.o.s., white circles mark the most significant substructures residing within R500, in the plane perpendicular to 
the l.o.s.. Spatial maps in columns (2)- (3)- (4) refer to the whole R500 region and, overplotted, is the ATHENA-like FoV for comparison. 
From top to bottom the haloes are ordered such that the value of /j, = av(< R500) / o'thermal (Eq.[6]) is decreasing. 
© ... RAS, MNRAS 000, [THIS] 



12 V. Biffi et al. 



vi.o.s. — kT phase diagram, color-coded by EM. The region 
of the clusters considered is always that enclosed within the 
projected along the line of sight through the simula- 

tion box. 

Despite the fact that all of these clusters show signif- 
icant non-thermal motions compared to the characteristic 
thermal velocity expected, they show different evolutionary 
states of such velocity patterns. 

Clearly, the maps show very different internal structures 
and the presence among these velocity outliers of both sub- 
structured haloes and more regular roundly-shaped ones, 
like halo 36 especially. The latter, in partiuclar, shows very 
regular features, no significant substructures and relatively 
low cooling time (as will be shown in Section [53)) . suggest- 
ing that the non-thermal velocities measured are not due 
to merging or infalling subhaloes but rather characterize in- 
trinsically the ICM. Similar is the case of halo 27, which 
shows a very regular morphology and structure from all the 
maps in Fig. [12] and, however, is charaterized by both large 
cooling times and significant non-thermal velocities. 

Such examples indeed represent the problematic case 
where no self-bound substructures are visible in the emis- 
sion maps, being probably dissolved but not yet fully ther- 
malized, which causes high-velocity streams to persist in the 
ICM. Here the usual classification of the halo as a regular, 
not disturbed cluster and the assumption of thermal motions 
to dominate the gas pressure support might be misleading. 

Nonetheless, we also find more trivial cases, like halos 37 
and 28, that show significant substructures in the emission 
maps, reflected in their velocity and temperature structure. 
This suggests a non-relaxed dynamical configuration, where 
the subhaloes are still in the form of self-bound objects with 
an own dark matter potential and are most likely responsi- 
ble for driving high- velocity motions. In fact, these clusters 
deviate significantly from the best-fit Lx —T relation, mean- 
ing that the kinetic energy of the infalling substructures has 
not yet thermalized and the boost in luminosity is not yet 
followed by the increase of temperature. Similarly, the high 
velocity measured for halo 29 might be related to a dis- 
turbed state of the cluster, for which a subhalo is clearly 
present from the emission map. 

In an observational case of this kind, though, clearly dis- 
turbed haloes can be identified and prominent substructures 
can be explicitely excluded from the analysis. 



5.3 Cool-core clusters in the simulated sample 

In Fig. [13] we plot the central cooling time, normalized to 
Hubble time, as function of centr al electron density, ne,o, for 
all the cluster in the sample. As in'Pratt et al.' ('2009*) , we de- 
fine the central values as those calculated at 0.03i?5ooLJ. The 
data reported have been derived from the three-dimensional 



10 While for global properties as Lx or T the difference in the def- 
inition of R500 is not crucial, here we also calculate the A = 500 
overdensity with respect to the critical density of the Universe, in 
order to compare as consistently as possible against the observa- 
tional data. In fact, the definition is particularly relevant in order 
to evaluate the values of cooling time and density at exactly the 
same radius, being not mean values within a region but rather 
local values. 
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Figure 13. Cooling time versus central electron density, with 
color code as in Fig.[9] Asterisks refer to simulated data, while 
black diamonds represent the observational sample investigated 
in Pratt et al. (2009) (P09). The dotted line marks the threshold 



adopted by P09 to define cool core systems, 
4 X 10~^ cm~^. 



h(z) ^ ne,o > 



profiles of the haloes, calculated directly from the simula- 
tion. 

Opposite to the overcooling problem of standard nu- 
merical simulations, we show with this Figure that no real 
cool cores are present in our sample, as a consequence of 
the inclusion in the simulations of AGN feedback. In fact, 
this has the effect of preventing the gas in the hot phase to 
remain in the central region long enough time to cool and 
form stars. Comparing to observational data, we do not find 
haloes populating the region of very short cooling times and 
high central electron density. 

Here, we adopt the same color code as in Fig. [9] in or- 
der to identify the most prominent velocity outliers previ- 
ously discussed. We note that the five outliers singled out by 
the velocity diagnostics (Section |5.2|) comprise both haloes 
with very long cooling times and halos for which the cooling 
time is only some tenths of the Hubble time. For the latter, 
although not classified as strong cool cores (CCs), cooling 
might be nevertheless important in the central region. With 
respect to our analysis, this can also be the case of sub- 
structured haloes, whenever the approaching subhaloes have 
not merged yet, namely the first passage through the clus- 
ter central region has still to occur. In such circumstances, 
the cooling time calculated for the cluster core can also be 
relatively small, since the innermost region has still to be af- 
fected by the interaction, while the luminosity and velocity 
fields already are. 

Among these haloes with relatively short cooling time, 
we find two haloes (37 and 28) out of the five clusters with 
significant non-thermal motions, which reside in the region 
above the best-fit Lx — T relation. This suggests that there 
is no real tension between their position in the Lx — T 
plane and the observational evidence of CC clusters occu- 
pying that region of the relation. The velocity diagnostics 
proposed can actually probe different aspects of the cluster 
state in addition to the usual dinstinction between CCs and 
morphologically disturbed clusters. 
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6 CONCLUSIONS 

In this paper we have presented the study of the ICM ve- 
locity structure and its impact on the Lx — T scahng re- 
lation and cluster status, from hydrodynamical simulations 
and synthetic X-ray observations of a set of galaxy clus- 
ters. Numerical simulations have been performed by means 
of the treePM/SPH, parallel code P-GADGET3, includ- 
ing a number of physical processes describing the bary- 
onic component with a level of detail never reached so far, 
namely cooling, star format ion and supernova-driven winds 
(jSpringel &: Hernauistll2003l ). chemical enrichment from stel- 
lar p opulation, AGB stars and SNe (iTornatore et alj|2004 
I2OO7I ). low-viscosity SPH (IPolag et all l2005l l. black-hole 
growth and feedback from AGN. The set of 43 simu- 
lated cluster-like haloes has been selected from a medium- 
resolution cosmological box by requiring M500 > 3 x 
10^^h~^ M0. The mock X-ray observations have b een ob- 
taine d with the virtual photon simulator PHOX (Biff i et al] 
12012) for both Chandra ACIS-I3 and the XMS spectrome- 
ter originally designed for ATHENA, and used here as pro- 
totype for any next generation X-ray spectrometer capable 
of reaching very high energy resolution. 

From the direct analysis of the gas component in the 
central part of the simulated haloes, in the projection plane 
perpendicular to the chosen line of sight, we find tight rela- 
tions among different definitions and selection criteria used 
to calculate the line-of-sight component of the gas velocity 
dispersion. In order to consistently compare and use results 
from simulations and mock observations, we always refer to 
cylinder-like regions which extend throughout the z coor- 
dinate of the simulated box, i.e. along the l.o.s. direction. 
In particular, we find that mass- and emission-measure- 
(EM-) weighted velocity dispersion calculated for the gas 
particles with i?5oo both correlate closely to the halo M500 
and, therefore, can be equivalent ly employed. Between the 
two definitions, we concentrate on the EM- weighted veloc- 
ity dispersion since it is more sensitively related to the gas 
X-ray emission. 

We also find that (t^qq^em probes fairly well the velocity 
dispersion in inner, smaller regions of the clusters, such as 
the region covered by a small FoV like that of ATHENA. As 
a fair approximation, we can safely assume the results to be 
representative of the whole region out to i?500 , even though 
single-pointing observations with a telescope of similar FoV 
of clusters at the considered redshift would only cover their 
innermost part. On the other hand, small-scale spatial fea- 
tures of the ICM velocity field could be possibly unveiled by 
the detailed comparison between the large-scale amplitude 
of the l.o.s. velocity dispersion and the value detectable with 
ATHENA-like observations. 

By applying our virtual telescope, PHOX, to the simu- 
lated haloes, the l.o.s. velocity dispersion is derived from the 
velocity broadening of the Ka iron complex around 6.7 keV 
in the synthetic X-ray spectra, provided the high energy 
resolution expected to be reached by such a spectrometer. 

Moreover, predictions on Gv provided by the simula- 
tion and detections mimicked with PHOX show a level of 
deviation of less than 20%, for 74% of the clusters in the 
sample. The largest deviations of observed velocities from 
expected values are significantly dependent on the compli- 
cated distribution of the gas velocities and thermal structure 



(see Fig.[7|). In principle, a more detailed modelling of the 
shape of the spectral emission lines, such as the iron com- 
plex used here, rather than the standard Gaussian fit, should 
permit to derive more accurate measures of the l.o.s. veloc- 
ity dispersion of the gas, accounting for multiple (thermal) 
components in the velocity field. 

Given the ideal case of a similar X-ray spectrometer at 
high resolution, we base our further investigation on the sim- 
ulation results. Indeed, this is further suggested by the good 
agreement between mock data and simulation and therefore 
we utilise the g^qq^em obtained directly from the simulation 
to constrain the impact of velocity structure onto X-ray ob- 
servables, such as luminosity and temperature. 

In principle, as soon as high-precision X-ray spec- 
troscopy will become available for real observations, we will 
be able to safely employ the detected values of the ICM 
non-thermal velocity dispersion for studies of real galaxy 
clusters as well. Possibly, this will be already explored, for 
instance, in the case of the brightest clusters by means of 
the upcoming ASTRO-H satellite, due to be launched in 
2014, for which a energy resolution of 7eV at 0.3 — 12 keV 
is expected. 

From the synthetic Chandra spectra obtained with 
PHOX, we were able to estimate observed bolometric lu- 
minosity, extrapolated to the entire X-ray band, and tem- 
perature, for the ICM residing within i^soo in the projection 
plane. The Lx — T scaling relation constructed from these 
mock data g enerally agrees with real observations, e.g. from 
IPratt et all JioOQiV By using the BCES linear regression fit- 
ting method (jAkrit as &: Bershadvll 19961 ) , we assume a linear 
relation between Lx and T in the log-log space and deter- 
mine the best-fit slope and normalization of the Lx — T for 
our 43 simulated clusters, find ing a shallower slo pe with re- 
spect to the results reported in IPratt et al.) (|2009l ). In the fit- 
ting procedure, we explore both the BCES (L|T) and BCES 
Orthogonal methods, where the former treats the luminosity 
as dependent variable and minimizes its residuals, while the 
latter implies both variables, L and T, as independent vari- 
ables and minimizes the orthogonal distances to the linear 
relation. 

As a step forwards, we include the information obtained 
about the ICM velocity dispersion along the l.o.s. in order to 
investigate the eflFects on slope, normalization and intrinsic 
scatter of the best-fit Lx — T relation. 

Interestingly, the exclusion from the original sample of 
the haloes with largest velocity dispersion (normalized to 
the characteristic thermal value associated to their temper- 
ature) has the effect of reducing the scatter in the Lx — T re- 
lation, as shown by the trend in Fig. [11] (bottom panel). The 
increasing trend of slope and scatter, in particular, evokes 
a dependence on the contamination of the sample by haloes 
which show significant fraction of non-thermal velocity, al- 
though the results would need a larger statistics to be more 
strongly confirmed, especially at the low-temperature end 
of the relation, which is not probed by our sample. 

Nevertheless, an interesting indication of our velocity 
diagnostics is that even haloes with fairly regular appear- 
ance or relatively low central cooling time can be charac- 
terised by high non-thermal velocities of the gas along the 
line of sight, deviating therefore from the expected scaling 
relations among global properties (e.g. the Lx — T scaling 
law) . Additional information on the velocity structure of the 
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clusters might therefore be helpful to constrain better their 
state. Observationally, for the clusters where the velocity 
measurements are indeed achievable, this can be used as 
complementary characterisation to the usual morphological 
or cooling-time-based approach. 

As a promising, future perspective, we expect high- 
precision X-ray spectroscopy to provide valuable informa- 
tion on the non-thermal velocity structure of the ICM in 
galaxy clusters along the line of sight, which can be safely 
assumed to trace the intrinsic gas motions of the cluster de- 
spite the effects due to projection and instrumental response. 
This will be certainly fundamental to correctly determine 
the total mass in clusters and characterise in more detail 
the deviation of real clusters from the hydrostatic equilib- 
rium at the base of X-ray mass estimates as well as the 
divergence from the self-similarity of the haloes expected 
from theoretical models. 
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APPENDIX A: TREATMENT OF CHEMICAL 
ABUNDANCES IN PHOX 

Dealing with galaxy clusters, the main emission components 
in the X-ray regime are constituted by bremsstrahlung con- 
tinuum and emission lines from heavy elements. 

Usually, simulations of the X-ray emission from hydro- 
dynamically simulated clusters assume the average value of 
one third the solar metallicity, consistently to the average 
metallicity measured for real clusters. With respect to previ- 
ous hydro-simulations, however, the run analysed here pro- 
vides full information on the chemical enrichment of the gas, 
tracing its composition with 10 elements, namely: He, C, Ca, 
O, N, Ne, Mg, S, Si, Fe. An additional field for each gas par- 
ticle accounts also for the total mass in all the remaining 
metals. 

Therefore, the available abundances of the individual 
elements. He included, are read and used to calculate the 
contribution of each element to the final spectrum. 

This is done, in PHOX, via the external package 
XSPEC, assumin g a single-temperature VAPEC (or, if de- 



sired, VMEKAlFl lMewe et al.1 1 19851 : iKaastra k Mewd Il993l : 
iLiedahl et al.fll995l ) model for each element, where the 



cor- 



responding abundance is fixed at 1 Z©. The total spectrum 
(S) of a gas particle, at its temperature T, will then be ob- 
tained as the sum of the continuum, due to H, plus the 
contribution of each element, weighted by its abundance 
in solar units (set, here, according to the values given by 
lAnders k Grevessd ([l989)). Explicitly: 

S{T, Ztot) = S{T, H) + S.Z. [S{T, Z.) - S{T, H)] (Al) 

The abundances of the VAPEC (VMEKAL) elements which are 
not explicitly treated in the hydro-simulation are estimated 
from the remaining mass in metals associated to the par- 
ticle, redistributing them according to the solar abundance 
pattern. 



APPENDIX B: THE SIMULATED SAMPLE 

In order to have a visual presentation of the sample discussed 
in this work, we show here emission measure (EM) maps for 
the 43 clusters, for the cylinder-like region along the line of 
sight, encompassed by i^soo around each cluster. 

This paper has been typeset from a file prepared 

by the author. 



See |http: / /heasarc. nasa.gov/xanadu/xspec/manual/XSmodelMekal. html | 
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Figure Bl. Visualization of the 43 most massive galaxy clusters in the simulated sample. The panels show the EM maps in the pL 
perpendicular to the line of sight, within the i^soo region of each cluster. 
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Figure B2. Continues from Fig. lBll 



© ... RAS, MNRAS 000, [THIS] 



